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Abstract 



We analyze the change in viewing angle of the holographic image for the on-axis 
hologram by varying an incidence angle of a plane wave. The tilted plane wave plays 
a role in a modulated carrier similar to a carrier signal of the off-axis holography, 
which makes the diffractive wave propagating to new optical axis direction. In the 
Presnel hologram, new viewing zone of the reconstruction image is formed, but its 
image shape is deformed along new viewing direction. The Fourier hologram enables 
to retrieve three-dimensional image with other perspective by a tilted plane wave. 
The sampled hologram with pixel structure generates high-order diffraction beams 
propagating to various directions. Collection of diffracted beams to a definite image 
spot effectively increases a viewing zone angle of the reconstruction image. We con- 
duct the numerical analysis for the sampled hologram showing high-order diffraction 
beams with various viewing zones. 

OCIS codes: (090.2870) Holographic display; (090.1995) Digital holography; 
(090.1970) Diffractive optics 
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1. Introduction 

The hologram imphes three-dimensional information with phase and amplitude of the prop- 
agating wave [1,2]. This makes it possible to retrieve the wavefront of original propagating 
wave by illuminating the hologram with reconstructing wave. However, in digital holography, 
the holography has difficulty in the recording and reconstructing processes of the hologram 
image because of high spatial frequency of the fringe pattern [3-5]. The hologram is the 
interference fringe of the recording wave and the object wave, and therefore, for retrieving 
the hologram image with sufficient viewing angle the spatial frequency gets to the extent of 
the wavelength of the light field. 

In holographic display, the hologram is sampled with sampling frequency of the device 
with pixel structure to avoid aliasing effects [6,7]. Present modulator encoding hologram 
fringe has much lower resolution than the required specification. Now, the limitation of 
the spatial frequency due to the device performance obstructs widening view angle of the 
hologram image, although the holographic display can reconstruct the image to the display 
size. For the realization of holographic display, solving the problem of widening viewing 
angle is essential. 

Current researches for improving the viewing zone have been carried out through the 
spatial multiplexing of the spatial light modulators [8-10]. Curved array of the devices 
effectively enhances the viewing angle of the reconstructed image, where from the stereogram 
representing each view, the large space-bandwidth product is still required. Since only the 
high-definition resolution of the digital image is sufficient, the massive hologram data for 
displaying holographic high-definition image is considered to be absurd. This also hinders 
the process of the digitized hologram data. 

In this research, we note that in principle, even the hologram with small space bandwidth 
has entire information of the object, and thus trying to find how to enlarge the viewing 
zone of the reconstructed image optically. First, for the Fresnel and Fourier holograms we 
investigate the possibility of the viewing angle change by varying an incidence angle of the 
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illuminating plane wave. The plane wave illuminating the hologram in other direction than 
that of the recording process generates diffracted beams emanating to the shifted position 
of the reconstructed image, which can make the viewing zone with another perspective of 
the image. Next, in the sampled hologram with pixel structure the method of extending the 
viewing zone angle by collecting high-order diffracted beams is proposed. Finally, we carry 
out the numerical analysis for the sampled hologram showing high-order diffraction beams 
with various viewing zones. 

2. Viewing angle analysis of hologram image 

A. Reconstruction process by a tilted plane wave 

The reconstruction process of the objective wave in the on-axis hologram illuminating by 
a tilted plane wave, exp(jA;oxsin^) is expressed as the product of a hologram component, 
a{x, y) exp[j0(a;, y)] and the illuminating wave. 



When the plane wave with unit amplitude is used in recording the on-axis hologram, a{x, y) 
and y) are real amplitude and phase of the object wave, respectively. Equation (1) can 
be rewritten as another form. 



This is a kind of the reconstruction form for the off-axis hologram made by the recording 
plane wave, ex^{jkQXsniO), where the coaxial reconstruction wave is assumed. Equation 
(2) can be also rewritten by varying together incidence angle values of the recording and 
reconstructing plane waves in the off-axis hologram. In a view of the carrier frequency 
holography, the on-axis hologram fringe is unique because it has a constant form as the 
modulating signal irrespective of the spatial carrier. Therefore, in two cases of the on-axis 
holography by an inclined plane wave and the off-axis holography, there is only the variation 
of the carrier frequency. 



a(x, y) exp[7'0(x, y)] exp{jkox sin 0) 



(1) 



a{x, y) exp {j [0(a;, y) + kgx sin 9]} 



(2) 
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Figure 1 depicts the reconstruction by illuminating the on-axis hologram with the plane 
wave at an incidence angle, 9. Here, we notice only the real image from the complex 
amplitude modulation in order to clarify our analysis. The viewing angle of the image by 
the coaxial plane wave is determined by means of the spatial frequency of the fringe pattern. 
Since the local spatial frequency is given by ^ fti ^^y^, the viewing angle, VL is double the 
size of ui. The higher spatial frequency leads the wider viewing angle of the image. In the 
tilted plane wave, since the exponential term of a carrier frequency generates the deflected 
beam of the retrieved image, the viewing zone of the object image changes to new ^new 
region in Fig. 1. 

We find that the tilted plane wave plays a role in a modulated carrier similar to a carrier 
signal of the off-axis holography, which makes the diffractive wave propagating to new optical 
axis direction. In the sampled hologram, the viewing zone of the reconstructed image has a 
limitation because it has the bandlimit signal to the Nyquist frequency to avoid the aliasing 
effect [11]. But, the viewing zone variation of the object image by a tilted plane wave is 
irrespective of the aliasing in principle. This case is different from the aliasing arising in 
sampling the modulated carrier signal to the hologram. Here, the modulating hologram 
signal of the object is already recorded in the sampled hologram, and only the inclined 
reconstruction wave related to the carrier signal is illuminated. Thus, by controlling the 
incidence angle of the reconstruction wave for the on-axis sampled hologram, we can get the 
reconstruction effect of the off-axis holography making a new viewing zone. 

B. Hologram image deformation by a tilted plane wave 

Let us investigate the change in the reconstructed image of the hologram transparency, 
g{x,y) by a tilted plane wave. The complex field of the image in the Fresnel propagation 
is developed as the convolution of the input field, exp[_7A;o(3:^sin^ -|- y sin (f)]g{x,y) and the 
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spatial impulse response, h{x, y; Zq) = exp{-jkoZo) exp \-j^{x^ + y^) 



2zo 



U(x,y;zo) = -r — exp 



-jTT-i^ +y ) 
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X exp 



-jJ^(^a:" + y" 
■^2zo^ ^ 



J J exjp[jko{x' sin 9 + y' sin (p)]g{x', y') 

dx'dy' (4) 
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The integral part can be calculated from the convolution operation of two Fourier transforms 
of the product term of the hologram and the plane wave and the quadratic phase term. 

1 
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[x + Zq sin 9,y + Zq sin Lp) 



*H{fxjY) 



(5) 



where H{fx,fY) is the transfer function in the Presnel diffraction. 

/ [ Gif'x + fxc, + fYc)H{fx - fx, fy - fY)dfxdf{. 



(6) 



By inserting Eq. (6) into Eq. (4), above equation is the angular spectrum representation 
of the field propagation in the Fresnel diffraction. If we substitute the variables, fx + fxc 
by fx and fy + fvc by fy and rearrange the equation, it becomes the expression for the 
complex field of the object shifted to the transversal axis. 



U{x + Zq sin 9,y + zq sin (p; zq) 



(7) 



The tilted plane wave causes the image reconstructed at the shifted position to the 
transversal axis. This analysis is applicable to the arbitrary incidence angle of the tilted 
plane wave, because the diffraction beam to retrieve an image can be regarded as the paraxial 
propagation along the new diffraction optical-axis. Figure 2(a) shows two-dimensional plane 
image vertical to optical axis retrieved at the shifted position by an inclined plane wave to 
X-axis. We find that in a plane object, new viewing angle of the image is generated, because 
only the lateral position of the reconstructed image is shifted without distortion of the object 
shape. 

However, in a three-dimensional object, a change of an image shape occurs. The three- 
dimensional object image can be calculated at various z values. 
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^U{x + zi sin 9,y + zi sin (p; zi) 



(8) 



The amount of shifted value to transversal axis depends on the z value, which invokes an 
image deformation along the depth direction. The retrieved image of the hologram composed 
of the point sources in x-z plane is illustrated in Fig. 2(b). The point image is obtained 
from the simple convolution expression about one-point hologram [12]. 



ex^{jkQX sin 0) exp 



J7^{x +y ) 

ZZq 



* h{x,y;zo) 



jj exp I J — [{x + zq sin 9)x' + yy'] | dx'dy' 5{x + zq sin 6*, y) 



(9) 



For convenience, the coefficient of each term shall drop. The first line of above 
equation is also presented as the reconstruction of the off- axis hologram as like 
exp |j —-^{x^ + y'^) + k^x sin 9 | y; Zq). The real image appears as the delta function. 
If sin 9 is written by the point image moves to Xp amount at x axis. Here, defiection beam 
angle does not exactly coincide with an incidence beam angle. Four real images are recon- 
structed at following locations, where the interval between two points in the axial direction 
is put to be c. 



(a; - I + Xpi, y]ZQ-^^ +5(^+^ + Xpi,y-zo-^ 

(c c\ ( c c 

^~ 2^ ^p2, y;zo + -j + 5 ix + - + Xp2, y;zo + - 



(10) 



where Xpi and Xp2 are the values, {zo — c/2) sin^ and {zo + c/2) sin^. 

As shown in Fig. 2(b), shearing transformation arises along x-axis. The image deforma- 
tion by a tilted reconstruction plane wave, exp[jko{xsm9 + y sin ip)] is described in terms of 
Affine transformation. 
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zsm9 



zsirnp 



(11) 



We find that the reconstruction plane wave incident on at a different angle from that of 
the recording process definitely induces the deformation of a retrieved image, which holds 
for both on-axis and off-axis holograms. Although the plane object image makes a new 
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viewing zone in terms of an inclined plane wave, three-dimensional image does not form 
new entire view in the Presnel hologram. That is, at a viewing direction observer sees the 
projection scene of reconstruction image by coaxial plane wave with respect to vertical plane 
of deflection axis, but deformed image shape. 

C. Viewing angle change for the Fourier hologram image 

Figure 3 illustrates the reconstruction process for the Fourier hologram by a tilted plane 
wave. The hologram plane is placed in the front focal plane of the convex lens of focal 
length, /. The field distribution of the image in the vicinity of the back focal plane, d is 
obtained from below equation [13,14]. 



where Az=f-d indicates the image depth distribution. From the similar interpretation of 
Eqs. (4)- (7) in the Fresnel hologram, we can get the expression for the complex field of the 
object image shifted to the transversal axis. 



Here, the shift quantity to the transversal axis has a constant value as fsm6 or fsinip 
irrespective of the image depth, and thus, as different from the Presnel hologram, the image 
deformation does not occur along the viewing direction. 

The tilted plane wave does not change the shape of the reconstructed image and only 
moves the image to lateral position, as depicted in Fig. 3. Therefore, the inchned plane 
wave can vary effectively the viewing angle of some object image reconstructed optically 
as a function of the incidence angle in the Fourier hologram. We clearly confirm that the 
perspective of the object image changes according to the propagation direction of diffracted 
beams. Considering the virtual image formation with opposite view of the real image, it 




(12) 



U{x + f sm9,y + f sin cp; Az) 



(13) 
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is apparent that even the hologram with extremely low spatial-frequency shall retrieve the 
whole view of the object image optically in a transparent object. This is possible because 
the hologram has entire information of three-dimensional object. 

We assume that the viewing angle of the object is widened by collecting above diffraction 
beams propagating to various directions. Figure 4 is a typical example of obtaining the wide 
viewing angle of the image by varying an incidence angle of the plane wave. Here, the 
deflection device should sufficiently collect diffracted beams to a definite image spot. If the 
device is far apart from the hologram plane, it will be a single type with respect to diffracted 
beams in Fig. 4. But if not, the device has a separate specification according to the beams. 
And for the large incidence angle, the correction of the optical path length is required as an 
amount, 2do(^-a — 1). 

' " ^ cos ' 



3. Expansion of viewing angle in sampled hologram 

The sampled hologram in holographic display generates high-order diffracted beams due to 
pixel structure. The diffracted beams propagate in the direction at an angle with respect 
to the optical ^-axis. We investigate the possibility of widening viewing angle by collecting 
diffracted beams to a definite image spot. Figure 5 illustrates the schematic diagram of 
forming the holographic image on the basis of our deduction. The optical element between 
the hologram component plane involving the Fourier lens and the image plane plays a role 
in collecting diffracted beams. That is, this device defiects diffracted beams to the image 
spot. 

We assume the sampled Fourier hologram, gs{^^ t\) by rectangular pixels in Fig. 6, which 
has the pixel interval and the pixel size Ap^ in the ^-direction, and p^, Ap^ in the 77- 
direction. 



00 00 



nj=— 00 nn=—oo 



(14) 



Ap^ Apr, ^ 

where rect() is a rectangular function. The amplitude distribution in front of the optical 
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element is represented by the convolution notation. 



Up{x, y) = [gs{x, y) * h{x, y; /)] exp 



* h{x,y;zi) 



(15) 



where zi is the distance between the hologram component plane and the deflection plane. 
We consider the normally incident plane wave with unit amplitude. The complex field in 
the image plane is obtained using above equations. 



U {u, v) = [Up{u, v)tp{u, v)] * h{u, v; Z2) 
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dxdy (16) 



The tp{x,y) is the amplitude transmittance of the deflection element, and the Z2 is the 
distance between the deflection plane and the image plane. 

The field distribution, Up{x,y) is simplified to the Fourier spectrum of the sampled 
hologram in case of the Zi equal to the focal length of the lens [13]. At first, for analysis of 
the property of diffraction beams we evaluate the Fourier transform of Eq. (14). 
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If the g{n^p^, n^Pr)) is a constant number, the summation part becomes the Dirac comb func- 
tion, which makes just the Fraunhofer diffraction of the rectangular pixel array. Inserting the 
value of the sampled hologram, this term turns into the Fourier spectrum. The summation 
term is rewritten as the Fourier transform through the Poisson summation formula. 
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From this, Eq. (17) describes the modulation of the periodic Fourier spectrum by the 
envelope of the sine function along the x- and y-axis. We briefiy depict the modulated curve 
by the sine function in Fig. 7. If the pixel interval p^^^ is put to be approximately equal to the 



9 



pixel size Ap^^^, the zeroth-order diffraction beam well implies the zeroth-order spectrum 
image. But the high-order images are centered at the minimum region of the envelope 
because the main lobe of the sine function has twice the width of other beams. This can 
be overcome by the phase-shift technology, which appropriately controls the position of the 
reconstructed image to the transverse axis. 

The complex field of the object in the deflection plane is distributed within each diffracted 
beam at the interval of y or Equation (18) can be expressed as another form as follows. 



E E i ^(C, ??) exp 

o=— oo /3=— oo V 



Pi Pv 



a=-oo /3=-oo 



(19) 



where the diffraction relation, sin 9a — aX and Pr, sin 9p — f3X, are applied. As stated 
in Section 2, above equation presents the diffraction of the g{C:V) object with the spatial 
frequency by illuminating with the inclined plane wave at an incidence angle 9a,i3 ■ There- 
fore, we find that high-order diffraction beams create the new viewing zone with different 
perspective of the object. 

Figure 8 is a system for collecting high-order diffracted beams to enlarge viewing angle. 
For convenience, only the first-order beam along the x-axis is displayed. The prism as an 
optical defiection device is used to bend diffracted waves to the same place of the image plane. 
The phase transformation of the deflection device is generally written by exp(— jA;a;sin6'). 
High-order diffraction beams should be deflected at the different angle according to the 
diffraction order. The amplitude transmission of this device has a separated value distributed 
with the space. Aw of the beams. 



tpix,y) = E^^P 
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+ ^exp 
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(20) 



The Wx^y is the shifted value equal to the beam space. Figure 9(a) exhibits the design of the 
deflection elements. 
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Prom the second term of the convolution in Eq. (16), we evaluate the Fourier transform 
of the amplitude distribution of the deflection device. 

22 



(21) 



2_^^d{u + Z2 sm ) * |smc(— — ) exp [^'' 2)}} 

+ 2_. \ o{v + Z2smeny) * |smc( ) exp [j^^ V^^^ 2;J J 

where the constant coefficient is omitted. In real system, since the space Aw of the diffraction 
beam in the deflection plane is very large, the sine function can be regarded as the delta 
function. Furthermore, considering the spatial distribution of deflection devices with an 
interval of w value, the function together with exponential term is expressed as like 5{u — 
Wx{nx + 1/2)). Equation (21) is reduced to be the following delta function. 

Y,s(u + Z2 sin - Wx + + + ^2 sin 6'„^ - Wy (uy + (22) 

Tlx 

The first term of the convolution in Eq. (16) is also rewritten as Up{u,v) * h{u,v;z2), 
where similar deduction in Section 2 is applied to obtaining the complex field, Ui{u, v) in the 
image plane. From this, the wave field distribution in the image plane can be interpreted 
by simplified convolution expression as below. 

a=-oo^=-oo \ Pi. Pv / 

* I X! ^ + ^2 sin 9n, - Wx (fix + (^^ + ^ny - Wy (uy + ^)) [ (23) 

I Tlx f^y \ 

In our system in Fig. 8, each deflection element operates only the diffraction beam with 
corresponding order, and thus the convolution is carried out between two terms with same 
order. If we select an adequate deflection angle corresponding to the diffraction angle of the 
beam, diffraction beams can be collected at a particular image spot. In case of equal distance 
of Zi and Z2 values, the deflection angle, sin 0^^^^^ becomes double compared to the diffraction 
angle, sin^^^^, as depicted in Fig. 9(b). The transverse location component compensates the 
shift quantity of the delta function to match with two terms of the convolution operation. 
We can flnally obtain all diffracted beams gathering in the zeroth-order region in the image 
plane. 
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j:UaAu',v') (24) 

a,/3 

Here, the deflection operation is exerted at the shifted location of the x-y deflection plane 
in Fig. 9(b). Therefore, in a view of the deflection plane, above equation implies the 
reconstruction by diffracted beams at a different angle. 

oo oo 

{u{x,y)exp\jk{x sin Oa + y sin Op)]} (25) 

a=-oo /3=-oo 

The wide viewing angle of the holographic image is realized only by collecting high-order 
diffraction beams with different perspectives. 



4. Numerical analysis of the sampled hologram 

The outline of the sampled hologram for the numerical analysis is illustrated in Fig. 10. 
The hologram is sampled in the rectangular pixel array, where the size of the pixel array 
is 200x200 and the interval of the pixels, p is set to be 4 //m. Based on the rectangular 
function in Eq. (14), each pixel is divided into several subpixels, which makes it possible 
to display high-order images in the entire area of the reconstruction picture. Here, the 
number of subpixels is put to be ten, and then, the total resolution of the device encoding 
the sampled hologram becomes to be 2000x2000. 

We prepared the Fresnel hologram with 200x200 resolution about the object of 'HOLO' 
letter in Fig. 11. The hologram was synthesized at z — 0.2 m apart from the object plane 
by using the recording plane wave with wavelength of 632.8 nm, where the border between 
pixels is one subpixel with 0.4 /im thickness and the pixel size, Ap becomes 3.6 /im. Each 
hologram data is recorded in the 200x200 corresponding pixel array. That is, the hologram 
data within a single pixel of the device has the same value. 

Figure 12 is the numerically reconstructed image by using above sampled hologram. For 
convenience, the complex modulation is considered. We flnd that high-order reconstructed 
images generate, whose intensity varies according to the envelope of the sine function. The 
ten high-order images up to fifth order are displayed in figure. The lower graph of Fig. 12(a) 
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shows the central cross-section spectrum along the x-axis. As discussed in Fig. 7, there is a 
misalignment between high-order images and the envelope of the sine function. Figure 12(b) 
is the reconstructed image for the hologram made by adding the phase-shift component. The 
zeroth order image is placed in the half main lobe of sine function, and the mismatch of 
high-order images is mostly improved. Here, half spatial bandwidth of the zeroth order 
diffraction is available. The board portion between pixels also causes the mismatch of the 
images. The envelope of the sine function and the image spectrum are placed with period 
of and y , respectively. The ratio of the pixel size, Ap to the pixel interval, p determines 
the number of the image in the one lobe of the sine function. We confirmed that the number 
of the image increases with decreasing the ratio. In above pixel structure, since the ratio is 
90%, there appears the mismatch of up to 10% per one order image. For well matching, the 
pixel size should be designed very close to the pixel interval. 

In Fig. 12, high-order images have the same plane picture, which is consistent with 
interpretation of Section 2. That is, diffraction beams only move the images to the trans- 
verse axis without their deformation. In order to investigate the object view of high-order 
images, the object image should be numerically reconstructed in the direction of the op- 
tical diffraction axis. Figure 13(a) illustrates the coordinate transformation to the optical 
diffraction direction. The each order image is reconstructed at the new Xo-yo coordinate. 
The images reconstructed in the several diffraction directions are appeared in Figure 13(b). 
The tilted angle is determined by diffraction order, where for visualization, large diffraction 
angle was adopted. The numerical simulation is carried out for each diffraction order in- 
dependently, where the interpolation operation is applied [15-17]. The reconstructed image 
is shrunk in the direction. The higher diffraction angle leads to the larger image 

contraction. The contraction of the images is natural because the reconstructed images can 
be regarded as the projection images to the tilted plane. Therefore, we find that especially 
in the Fourier hologram the object view varies according to the diffraction order, as depicted 
in Fig. 3. The variation of the viewing angle shall be clarified by the numerical simulation 
of three-dimensional object. 
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Finally, we can extend the viewing zone angle of the reconstructed image by collecting 
diffraction beams to the main lobe of the sine function. For this method, all of diffraction 
beams should have the same intensity to form the uniformly extended viewing-zone. The 
optical attenuator is one of methods for adjusting each beam intensity and correcting the 
optical path length. Furthermore, the spectrum intensity of a particular order can be ar- 
bitrarily controlled by the design of the diffracting grating. Therefore, the realization of 
this holographic display depends on the design of the pixel structure sampling the hologram 
data. 

5. Conculsion 

The illuminating plane wave to the Fresnel hologram at a different angle from the recording 
process leads to the formation of new viewing zone accompanied by the image deforma- 
tion along viewing direction. The Fourier hologram makes it possible to retrieve three- 
dimensional image with other perspective by a tilted plane wave. In the sampled hologram 
with pixel structure, the viewing zone angle of the reconstructed image is effectively in- 
creased only by collecting high-order diffraction beams with different viewing direction. The 
design of the pixel structure creating high-order diffraction beams uniformly distributed is 
important for this holographic display. Further study for making various types of diffrac- 
tion beams is required. This becomes one of the methods for overcoming the low spatial 
frequency making the narrow viewing angle. 
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Fig. 2. (a) Viewing angle change of the plane object image and (b) deformation of the image 
along depth direction reconstructed by a tilted plane wave incident on the on-axis Presnel hologram. 
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Fig. 4. Typical diagram of obtaining a wide viewing angle of the holographic image by varying 
an incidence angle of the plane wave. 
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Fig. 5. Schematic diagram of describing the collection of diffraction beams emanated from the 
sampled hologram. 
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Fig. 6. The pixel structure of the sampled hologram. 




3lz 2Az Az ^ ^ 23Lsr Site 

Vx Vx Vx Vx Vx Vx 

Ar, Ap, Ap^ Ap^ Ap^ Ap;, 



Fig. 7. Modulation of the periodic Fourier spectrum by the envelope of the sine function along 
the X-axis. Right inset is the diffraction beam images. Left inset is the resized figure to the 
logarithmic scale. 
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Fig. 8. Illustration of holographic display system for collecting high-order diffraction beams to 
the image plane by using the optical deflection device. 




Fig. 9. (a) Alignment of the optical deflection devices and (b) the deflection of diffraction beam. 
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Fig. 10. Outline of the sampled hologram for the numerical analysis. 




NOLO 



(a) 



(b) 



Fig. 11. (a) Sampled hologram with 200x200 resolution prepared from (b) the object of 'HOLO' 
letter. 
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Fig. 12. Numerically reconstructed images and image spectrums for (a) the sampled on-axis 
hologram and (b) the hologram made by using the phase-shift technology. Central cross-section 
spectrum along the a;-axis is presented. 
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Fig. 13. (a) Geometry of the coordinate transformation to the diffraction direction, (b) The 
several object views reconstructed at the tilted axis. 
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